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ABSTRACT 

We investigate the extent to which the pure magnification effect of gravitational lensing 
can be extracted from galaxy clustering statistics, by a nulling method which aims to 
eliminate terms arising from the intrinsic clustering of galaxies. The aim is to leave 
statistics which are free from the uncertainties of galaxy bias. We find that nulling 
can be done effectively, leaving data which are relatively insensitive to uncertainties 
in galaxy bias and its evolution, leading to cosmological parameter estimation which 
is effectively unbiased. This advantage comes at the expense of increased statistical 
errors, which are in some cases large, but it offers a robust alternative analysis method 
to cosmic shear for cosmological imaging surveys designed for weak lensing studies, or 
to full modelling of the clustering signal including magnification effects. 

Key words: cosmology: cosmological parameters - gravitational lenses - large-scale 
structure of Universe 



1 INTRODUCTION 

Cosmic shear has long been recognised as a potentially very 
powerful tool for determining the properties of the Universe. 
Its sensitivity to cosmological parameters and its clean de- 
pendence via gravity on the mass distribution make it an 
attractive option for determining, amongst other things, the 
properties of Dark Energy, and it also opens up the prospect 
of testing beyond-Einstein gravity models. For recent re- 
views and summaries, see Munshi et al. (2008); Hoekstra & 



Peacock et al. (2006); Albrecht et al. (2006) 



With the advent of very large surveys directed towards 
large-scale weak lensing effects, such as Pan-STARRSl, the 
Dark Energy Survey, and in the future, LSST, Euclid and 
WFIRST, focus has shifted towards the systematic errors 
rather than the raw statistical power of weak lensing. For 
cosmic shear, there are challenges related to the precise mea- 
surement of galaxy shapes (e.g. K itching, Taylor Heavens 
lAmara Refre gier 2008 | [Voigt Bridle||2010| |Bern- 



2008 



stein 



2010| [Kitching et al.||2010 ), and physical effects such 



as intrinsic alignments (e.g. [Heavens, Refregier Heymans 



Metzler 2000; Crittenden et al. 2001; Cate- 
Jing 2002; Mackey^ 



Ian, Kamionkowski Blandford 2001 

White Kamionkowski||2002| [Hirata Seljak 2004; Bridle 



k Abdalla|2007| [Schneider k Bridle|2009| ) which can mimic 
to some degree the effects of weak lensing. These effects will 
increase the errors on recovered parameters, so it is timely 
to investigate magnification (or amplification) - another fea- 
ture of weak gravitational lensing, which, although not as 
statistically powerful as cosmic shear, provides independent 
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cosmological information. Given that the effects of chang- 
ing Dark Energy properties and changing the gravity law 
are rather subtle, it is highly desirable to have independent 
techniques which are not subject to the same systematic 
errors. 

Weak lensing induces small changes to object sizes, 
which, through the conservation of surface brightness, leads 
to flux changes which alter the number density of galaxies 
when a flux limited survey is considered. As a result of this 
amplification, in concert with the changes in positions, the 
clustering pattern of galaxies is altered. In regions where 
there is magnification (as opposed to de-magnification), the 
number density is increased or decreased depending on the 
slope of the number counts, as it depends on how many 
sources are promoted into the survey by the lensing, and' Van| 
Waerbeke^ (2010) has shown how magnification and shear 
can be used together to help cosmological parameter esti- 
mation. Recent studies (jScranton et al.|2005 Hildebrandt et 



al.|2009 Menard et al.|2010 ) have demonstrated the magni- 
fication effect with cosmological surveys, and ' Zhang Penj 
( 2005 2006 ) have shown how flux information may be used 



to isolate lensing magnification. 

The main complication which has to be dealt with is 
that there is a strong intrinsic clustering of galaxies, which 
contaminates the clustering signal due to weak lensing. In- 
deed, to call it a contaminant is a misnomer, as it is much 
stronger than the lensing signal. The intrinsic clustering of 
galaxies is not as easy to model theoretically as the lensing 
signal, as the former depends on galaxy formation efficiency 
(normally distilled into galaxy bias), whereas the lensing sig- 
nal depends only on the matter distribution. One approach is 
to model the entire clustering system, including lensing and 



bias, and this has been the approach of Bernstein (2009); 
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Yoo (2009); Joachimi k Bridle (2010). However, since the 



main advantage of lensing is its independence (at least in 
its dominant effect) on complex astrophysics which may be 
difficult to determine with precision, it is appealing to con- 
sider an alternate strategy, where one removes the intrinsic 
clustering part of the signal, leaving, ideally, the pure lens- 
ing signal. If this can be achieved, then confrontation with 
theory becomes much more straightforward and robust. The 
purpose of this paper is to investigate to what extent this 
can be achieved, using a 'nulling' technique similar to that 
which has been proposed to deal with intrinsic alignments 
in cosmic shear. Nulling can be achieved because we know 
(for a given cosmology) how the source of the magnification 
signal depends on redshift, and we can exploit this to find 
combinations of correlations which do not include intrinsic 
clustering terms. 

The problems of dealing with intrinsic alignments in 
cosmic shear and intrinsic clustering in cosmic magnifica- 
tion are in some respects similar, but in others different. In 
both cases the observable quantity (the ellipticity in the case 
of shear, the number density in the case of magnification) 
depends to first order on the addition of two components, 
so any quadratic quantities such as the correlation function 
or power spectrum of the observables have four terms. The 
main differences are that in the cosmic shear case, the con- 
taminants are relatively small, but poorly-known theoreti- 
cally, whereas in the magnification case, the contaminating 
terms are large, but in principle measurable. 

The structure of the paper is as follows. In section [2] 
we review the theory of cosmic magnification. In section 3 
we develop the nulling theory, and in section 4 we present 
results, with conclusions given in section 5. 



2 COSMIC MAGNIFICATION AND 
AMPLIFICATION 

The distortion of light bundles leads to a mapping from the 
source plane to the image plane which is described by a 
symmetric amplification matrix which is decomposed into a 
convergence k,, and two shear terms 71 and 72, which de- 
scribe distortions along the coordinate axes and at 45° re- 
spectively: 



The solid angle elements are also magnified, thus lowering 
the number density and resulting in the factor (Narayan 



A: 



^ + 71 
-72 



-72 



Size changes (and hence flux changes, by Liouville's theo- 
rem) are related to the inverse of the Jacobian of the trans- 
formation, which to linear order is 

This is the magnification or amplification. Let the unmagni- 
fied average number density at position r above a flux limit 
/ be n(> /, r), and we assume that near the flux limit of a 
survey it has a slope a, i.e. 

n(>/,r)<x/-° 

then the expected observed galaxy counts, neglecting intrin- 
sic clustering, for given magnification /i(r) will be 

^(> /,r) = n{> /,r)/i(r)'^~\ 



1989). Note that the slope may be dependent on r (only 
via r = |r| if the survey selection is isotropic within the 
observed area), and that the power- law need not extend over 
a wide range of /; it is the local slope at the flux limit 
which is required. If we include intrinsic clustering of the 
sources, characterised by a fractional overdensity 6g, then 
the observed number density of galaxies at position r is given 
by (dropping the > /) 

n(r) = n(r)Mrr-Ml + <5,(r)] 

where 6g is the intrinsic fractional overdensity of galaxies. 
If we define b as the bias parameter, here in real space, so 
6g = b6, where S is the fractional mass overdensity, then 
linearising gives 

n(r) = n(r) [1 + Xk + b6{r)] , 

where A = 2(a — l) may be a function of r, as may b, through 
time-evolution. A should be relatively easy to measure from 
sufficiently deep photometric catalogues, so we initially re- 
gard it as fixed and known. The convergence is 



dr Fxir, r') 5{y') 



where FK{r,r) = r{r — r)/[ra{r)] for a fiat Universe (for 
r > r , otherwise), and A = 3ifo^m/(2c^). We wifi assume 
fiatness, but this can be easily generalised. Thus since k is 
linear in 5, we have to linear order that the fractional number 
overdensity A = n/n — lis: 

A(r) = X{r)K{r) + b{r)S{r). 

Let us consider the two-point statistics of A, for a tomo- 
graphic survey divided into ntins shells in estimated redshift 
(normally coming from photometric redshift estimates). We 
will frame this discussion in real space, but will switch to 
Fourier space in section [3] The two-point cross-correlation 
function of the number overdensity in two shells i and j is 
(with obvious notation) 

(AiAj) — XiXj{KiKj) -\- bibj{SiSj) 

+ Xibj{KiSj) -\- biXj{SiKj). ij = l,...nbins 

If we wish to have a clean cosmological test, then it is the 
first term which is of most interest, as it does not contain 
galaxy bias, which is not precisely known. If we are to use 
it, then we need to remove or model the other terms. In this 
paper we choose the former strategy, which can be achieved 
by nulling out the remaining terms. This is very similar to 



what has been proposed by [Joachimi Schneider (2008) 
and J oachimi Schneider] ( |2009 t for removal of intrinsic 
alignment (I A) terms in cosmic shear, where the second term 
is analogous to the so-called II term, and the cross-terms 
are GI terms in cosmic shear (and only one survives in that 
case). We turn to how to do this in the next section. 

There are two main differences between magnification 
nulling and lA. One is the magnitude of the effect. GI is gen- 
erally subdominant to the desired GG shear signal, whereas 
here the terms far exceed the bias-independent term. 
See the top panels of Fig. 1. The other difference is that we 
have rather poor theoretical understanding of the GI term, 
whereas we can in principle make some measurements of 
bias of galaxies. 
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3 NULLING 

As with the shear nulUng method of Joachimi & Schneider 
(2008), we cross-correlate the overdensity of galaxies in a to- 
mographic bin labelled by i with a weighted average of the 
overdensities in all other tomographic bins. Writing this av- 
erage as D^{0) in the continuum limit, where is an angular 
position on the sky. 



D\e) 



dr W\r)/\{Y). 



If we choose W^{r) to be zero for r close to Vi (the centre 
of the tomographic bin i) , then the intrinsic clustering term 
(a product of h5 terms) will be small. This is analogous to 
the downweighting ( Heymans k, Heavens||2003| or removal 
of ( [King Schneider||2002[ |2003t close" pairs to remove the 
II term in cosmic shear. 

We also want to zero out the cross terms {k5 terms). 
One cross term is 

poo nr 

Ah{ri)5{Yi) I drW\r)\{r) I dr FKir,r)6{r), 
Jo Jo 

and we choose such that this is zero. Reversing the in- 
tegration order (and dividing by A), 







where 



b{ri)S{ri 



POO POC 

•I 



drW\r)X{r)^^^^^^6{r) 



ra{r') 



g\r) = I dr'W\r')\{r') ^^ /\ 



(1) 



(2) 



If we choose g^{r) = on scales where the density corre- 
lations are significant (i.e. at r :^ r^), then the expectation 
value of equation[l]will be nearly zero. With r^ at the median 
redshift of the bin, we choose 

g\r.) = 0. 

Since K,{ri) is influenced only by r < r^, we null the other 
cross term by simply requiring 



W\r) = 0; 



r < n 



There are many weight vectors which will satisfy this equa- 
tion, so for maximum statistical power, we follow poachimi| 
Schneider| (2009) and include all orthogonal weight vectors 
which satisfy the constraint. Indeed the weights are identi- 
cal to the intrinsic alignment case if the slope of the number 
counts does not vary with redshift, but is modified if it does. 
For the details of how this is done, the reader is directed to 



Joachimi & Schneider (2009). Note that the weights require 



the choice of a fiducial cosmology, and the weights must not 
be changed as a search of parameter space is undertaken. 
Clearly, if the choice is very poor, the nulling will be subop- 
timal, but there would be no difficulty in iterating the pro- 
cess, where the best-fit cosmology is used to define optimal 
nulling weights, and the process repeated. In this paper, we 
assume the correct cosmology for the weights, for simplicity. 

To compute the weighted cross-power spectrum, we 
note that we are cross-correlating two fields 



A(rj)«;(ri) + b{ri)5{ri) 



(3) 



f 



drWl{r) [A(r)At(r) + 6(r)^(r)] 



s > 0, 



where the first term averages the overdensity over the shell. 
In the second term, we allow for a variety of weighting func- 
tions s — 1, . . . Uhins — 1, as there are Uhins — 1 orthogonal 
weighting functions which satisfy the constraint equation [2] 
We can write these in the form 



f 

Jo 



where 
^o(^) 



dr ql{r)6{r) 



AX{r)^^^-^--^ + b{r)T{r - r^ 



nair) 



(4) 



qs{r) = Ag\r)—- + h{r)W\r) s>0 
a(r) 



where T is roughly a top-hat of unit area centred on the bin 
centre r^. This can be refined to reflect the averaging over 
the bin width with a varying number density. 



For such fields, the Limber approximation (e.g. Bartel- 
mann Schneider|200H section 2.84) gives the cross-power 
spectrum as 



f 

Jo 



where the last argument of the matter power spectrum P 
indicates its time-dependence. 



3.1 Fisher analysis 

The sensitivity to cosmological parameters comes from the 
matter power spectrum and its evolution with time, and the 
distance-redshift relation 



r{z) 



r dz' 

Jo H{z^) 



where H{z) = Ho [(1 - Q)(l + zf + Qde + ^m{l + zf] 
if the Dark Energy has present density parameter Qde- 
and Q are the present matter and total density contribu- 
tions, and we ignore radiation. Note that later we allow the 
Dark Energy density to evolve, so this equation is modified 
in the standard way. 

For this initial study, we compute the covariance of 
power spectrum estimates assuming gaussian statistics, but 
note that this would need to be estimated from numerical 
simulations (e.g. 'Kiessling, Heavens Taylor|2010| for real 
applications, and would include non-zero cross terms. The 
covariance of powers averaged in bins of width Ai is 



27T 



l\-^OO.l^ss' 



where AQ is the solid angle of the survey, and 
dr 



-^Os'£-^sO,i 



f 

Jo 



-ql{r)ql,{r)P{i/r;r) -\- shot noise 



is a function of the bins i and j. For i ^ t the covariance is 
zero if we assume all-sky coverage. For practical cases we use 
the usual scaling with the sky fraction /sky if, as is the case 
here, most of the signal comes from relatively small scales. 
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Finally, we compute the Fisher matrix (Tegmark, Tay- 
lor &: Heavens|1997| for the cosmological parameters 0^: 



£ bins 



(5) 



and the summation convention is assumed for a, b, which are 
indices which run over all the PqI modes. 

The Fisher matrix also allows us to compute analyti- 
cally the parameter biases arising from incorrect parametri- 
sation of the data, assuming a gaussian likelihood surface 
(see e.g. Knox, Scoccimarro Dodelson||1998 Taylor et al. 
|2007|: 



where 6ip^ is the error made in the fixed model parameter 
^/3, and Fe^p is a pseudo-Fisher matrix of the same form as 
equation [5] but with one replaced by i/j. 



4 RESULTS 

The situation which we investigate is as follows. The galaxy 
bias is scale-independent, but redshift-dependent: 

b(z) = l + mz + Ah 

with m = 0.1 and Ab = 0.2. For the analysis, we incorrectly 
model the bias as a constant, independent of z, and treat 
it as a parameter in the modelling, with a gaussian prior 
of width 0.1, centred on unity. We also have 6 cosmological 
parameters (scalar spectral index, density parameter in mat- 
ter and baryons. Dark Energy equation of state parameter 
'^0 = p/{pc^)i fluctuation amplitude and Hubble parame- 
ter, with fiducial values ris = 1, Qrn = 0.25, = 0.05, 
it;o = — 1, cTs = 0.8, h — 0.7. We take a number count slope 
a = 1.5, independent of redshift. Note that this is a sim- 
plification to illustrate the method; in practice a depends 
strongly on magnitude limits and redshift. See, for example, 
the discussion in Appendix C of ) Joachimi Bridle (2010). 

For the survey, we assume a Euclid-like basic survey 
design of 20,000 square degrees, but initially we will ignore 
photometric redshift errors for illustration. We assume 35 
galaxies per square arcminute, with a redshift distribution 
n{z) oc exp[— (z/zo)"^"^], and a median redshift of 0.9. The 
distribution is truncated at z = 3, and 20 redshift bins are 
used, chosen such that the number in each bin is constant, 
with 50 angular power spectra computed on a logarithmic 
grid between / = 10 and 3000. In all cases, the auto-power 
spectra are not considered, to suppress the 66 term. 

In Fig.[l]we show in the lower panels the effect of nulling 
on the cross-power spectra. We see that if we have excellent 
redshifts for the galaxies, then the intrinsic clustering term 
{K,i6s) can be nulled out very effectively, reducing its rela- 
tive amplitude by nearly two orders of magnitude so that it 
is well below the {huKs) term (here we use a compact no- 
tation, meaning the cross-power of the k component of the 
overdensity in bin i with the relevant component of one of 
the weightings s associated with that bin) . The nulling is not 
perfect, because of the finite width of the bins: the nulling 
works at the median redshift of the bin, but is not perfect 
away from the bin centre. We see that we can reduce the 



contamination to a level generally below the {huhis) term, 
but it is still non-negligible. This figure suggests that some 
modelling of the bias is still needed, however the accuracy 
required is reduced. 

4.1 Modelling vs nulling 

The two general approaches we can take are to model the 
data, including a bias model, or to attempt to null out the 
terms depending on bias, and to model the remaining terms, 
which depend only rather cleanly on cosmological parame- 
ters. Our results of the effectiveness of nulling indicate that 
the latter approach is too ambitious, as the residual hi6 sig- 
nal is not small enough to be ignored, and it must be mod- 
elled. Having said this, its sensitivity to bias is much smaller 
than without nulling, so there are some advantages to this 
approach. What we show in this section is that if the mod- 
elling of bias (and its evolution) is imperfect, then modelling 
the complete clustering signal can lead to large biases in 
cosmological parameter estimates (in terms of the statisti- 
cal error), whereas with nulling, the parameter biases are 
typically much less than the statistical errors. On the other 
hand, the statistical errors with nulling may be consider- 
ably larger than in the non- nulled case, so nulling gives us 
a rather conservative, but robust analysis. 

We perform various analyses. In Version 1 (nulled), we 
perform nulling to remove the k6 terms as far as possible. 
Here we assume that the nulling is perfect, and there is no 
signal left except for the pure term. Because of the finite 
width of the redshift bins, the nulling is not perfect, and 
leaves a residual Hi6 term, which leads to some bias in the 
parameter estimates, some of which remain unacceptably 
large, and very large statistical errors. Thus we conclude 
that even when the photometric redshifts are assumed to be 
known exactly, one cannot ignore the residual Hi6 contami- 
nation in the nulling, and it has to be modelled. 

In Version 2, we model the residual Hi6 terms. In the 
standard case, we are essentially modelling the entire signal, 
without attempting to reduce the k,6 term by nulling. We see 
from Table ^ that for the bias model we investigate, the in- 
correct assumption for the bias evolution (i.e. that there is 
none) leads to residual biases which are typically quite large 
compared to the statistical errors. The extent of the biases 
is dependent on how wrong our bias model assumptions are 
- if we had chosen the correct parametrisation for the bias 
evolution, we would not expect a bias in the recovered pa- 
rameters. However, it is not clear how well we will be able to 
constrain the bias. We can use the microwave background 
radiation and the measured galaxy power spectrum, if we 
assume a gravity model, or use higher-order statistics (e.g. 



Verde et al. (2002)), but it will probably be difficult to con- 
strain the bias very well. In the nulled case, we remove most 
of the K,6 signal, but we model the residual contamination. 
Since the contamination (which depends on the galaxy bias) 
is relatively small, the error in the assumed bias evolution 
has a much smaller effect, and the biases in the parameter 
estimation are negligible in comparison with the statistical 
errors. The main reason for this is that the statistical er- 
rors have grown substantially (by factors of 3 — 12), but 
we see from the table that for most parameters, the biases 
are in absolute terms smaller than in the no-nulling case. 
These results are illustrated in Fig. [2] for a few of the cos- 
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Figure 1. Nulling performance for an example set of power spectra for foreground bin i = 3. Magnification correlations {ntx terms) are 
given as solid curves, magnification-clustering correlations {i^S) as dotted lines. The top panels show from left to right the cross power 
spectra of bin 3 ~ 0.46) tomographic bins j = 4, 5, 6 (with central redshifts z = 0.53, 0.59 and 0.65. The bottom panels show nulled 
power spectra for the first three orthogonal weight vectors s = 1,2,3 (note that in the variant of nulling employed in this work no 
particular ordering is associated with s, which each represent a linear combinations of signals from different tomographic bins. See text 
for more details). 



Table 1. Resulting Icr marginalised statistical errors (cr) and parameter biases (6, not to be 
confused with galaxy bias) from Fisher matrix predictions for different analysis strategies. 
Version 1 corresponds to treating i^hz as the signal and ignoring the contamination by the 
term, which acts as a systematic error. Version 2 corresponds to treating the sum of tx5 and 
the imperfectly nulled term as the signal, marginalising over the assumed constant galaxy 
bias in the latter contribution. Here a systematic error arises because the bias model for the 
1^5 term is not correct, but we see that nulling very effectively reduces this systematic to a 
low level. 



parameter 


Version 1; nulled 


Version 2; standard 


Version 2; nulled 




a 


b 


cr 


b 


a 


b 




0.463 


0.138 


0.012 


0.000 


0.078 


-0.021 




0.516 


-0.222 


0.009 


0.020 


0.120 


0.038 


h 


19.650 


6.238 


0.134 


0.382 


0.554 


-0.059 


Us 


3.888 


-1.165 


0.051 


-0.165 


0.211 


0.033 


Qb 


2.324 


0.760 


0.011 


0.024 


0.030 


-0.006 


wo 


4.780 


-0.239 


0.088 


0.202 


0.267 


0.008 



mological parameters. The nulled statistical errors (orange, 
shaded regions) are much larger than modelling the full sig- 
nal (blue, solid regions), but the bias is typically smaller. 
Note, however, that the performance of nulling compared 
with full modelling depends on how wrong the assumed bias 
model is. If our model is poor, then nulling can fare better 
in absolute terms. We turn to this and related issues next. 



4.2 Sensitivity to assumptions 

In Fig. [3] we illustrate how sensitive the errors are to an as- 
sumed prior on the galaxy bias. The prior is assumed to be 
gaussian, centred on unity, with r.m.s. at, labelling the x 



axis. Each panel illustrates the errors on a cosmological pa- 
rameter, marginalised over the others. General features are 
that the nulled statistical errors are always larger than the 
fully-modelled case. In most cases, the nulled bias is smaller 
than the fully-modelled case (for the fully-modelled bias 
is relatively small, but this seems to be anomalous). The 
nulled bias is always negligible in comparison with the nulled 
statistical errors, but the opposite is the normal situation for 
the fully-modelled case. This is accentuated when the prior 
on the galaxy bias is very small - we are essentially enforc- 
ing an incorrect model for the bias. When we increase the 
width of the galaxy bias prior, it is still unable to model the 
bias evolution (as it assumes there is none), but at least it 
allows the assumed constant bias to increase to closer to the 
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Figure 2. la constraints from the and residual hiS signals, including marginalisation over 
the galaxy bias. Blue filled contours correspond to marginalisation only (3rd column of Table 
1); orange shaded contours correspond to marginalisation + nulling (4th column of Table 
1). Parameters not shown are marginalised over. The cross in each panel marks the fiducial 
parameter values. 



average bias of the galaxies (at redshift 1, the galaxy bias is 
1.3). 

In Fig.|4]we see the effect of varying the slope m of the 
galaxy bias evolution model. This illustrates a main conclu- 
sion of this study, which is that nulling is particularly ef- 
fective if one makes incorrect assumptions about the model 
parametrisation. Here m is fixed, and not marginalised over, 
and all the panels in the figure show that the worse the 
parametrisation, the greater the biases in the parameter es- 
timates. The solid grey (red) lines show very large param- 
eter biases if the assumed slope is significantly too low or 
too high, and in some cases the fully-modelled case fares 
extremely badly. Inspection of the lower right figure shows 
that the systematic bias on the Dark Energy equation of 
state can far exceed the statistical error in the nulled analy- 
sis. Interestingly, this panel shows that nulling is very effec- 
tive in this case, showing a negligible systematic bias, even 
of the galaxy bias evolution assumptions are very wrong. 
This insensitivity to the parameterisation of the galaxy bias 
evolution seems to be a fairly general feature of nulling, and 
is its major strength. 



In the case of negligible photometric redshift errors it 
is sufficient to remove redshift auto-correlations from the 
analysis in order to avoid contributions from the 6S term. 
For a finite photometric redshift uncertainty however, SS is 
non- vanishing also for cross-correlations between neighbour- 
ing redshift bins, and we include this effect in our study. 
To reduce this further potential source of parameter bias, 
we adopt the procedure proposed by Joachimi Schneider] 
( |2009 ) who downweight correlations between adjacent bins 
in the transformed power spectra with a Gaussian centred 
on the mean redshift of the foreground redshift bin, choos- 
ing a width that is determined by the photometric redshift 



scatter (see Joachimi Schneider ( 2009 ) for details) 



In Table |2] we see the effect of a more realistic photo- 
metric redshift error of 0.03(1 + z), which is the goal for 
Euclid. We still find that nulling introduces negligible cos- 
mological parameter biases, but the statistical errors are in- 
creased, considerably in some cases. Perhaps only for the 
Dark Energy equation of state parameter is there a strong 
case for preferring nulling over full modelling, for the galaxy 
evolution model considered here. 



Cosmic magnification: 



nulling the intrinsic clustering signal 7 




0.22 
0.2 
0.18 
0.16 
0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 



0.05 



a; fully modelled 
b; fully modelled 
a; nulled 
b; nulled 



0.1 



0.15 



0.2 




0.3 
0.25 

0.2 
0.15 

0.1 
0.05 





0.05 



a; fully modelled 
b; fully modelled 
a; nulled 
b; nulled 




0.1 



0.15 



0.2 



Figure 3. Statistical la errors (black lines) and biases (red/grey lines) for different cosmological parameters as a function of the prior 
on the galaxy bias parameter, keeping all other bias model parameters at their fiducial values. Solid curves correspond to the results for 
the fully-modelled power spectrum analysis including marginalisation over the galaxy bias, dotted curves to the combined nulling and 
galaxy bias marginalisation approach. 

Table 2. Like Table^ but for a Gaussian photometric redshift scatter with width 0.03(1 + 2;). 



parameter 


Version 1; nulled 


Version 2; standard 


Version 2 


nulled 




cr 


b 


cr 


b 


cr 


b 




0.426 


0.080 


0.012 


-0.003 


0.136 


-0.063 




0.468 


-0.057 


0.009 


0.018 


0.165 


0.097 


h 


16.788 


-0.986 
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-0.148 


Us 
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0.079 


0.058 


-0.202 
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0.019 




1.992 


-0.091 


0.034 


0.046 


0.694 


-0.042 


wo 


4.338 


0.246 


0.087 


0.211 


0.260 


0.023 
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Figure 4. Statistical la errors (black lines) and biases (grey/red lines) for different cosmological parameters as a function of the slope 
m of the galaxy bias model, keeping all other parameters at their fiducial values. Solid curves correspond to the results for standard 
analysis including marginalisation over the galaxy bias, dotted curves to the combined nulling and galaxy bias marginalisation approach. 



In Fig. [5] we see the effect of varying the zero-redshift 
galaxy bias, bo. The true 6o is 1 + A6, where Ah is the offset 
shown in the figure. Recall that the bias is assumed to be 
constant, with a gaussian prior of width 0.1 centred on bo = 
1, and the true bias is 1 + 0.1^ + Ab, so the average bias for 
galaxies is larger than 1 + Ab. Hence for Ab > 0, the model 
finds it hard to match the average galaxy bias for the lensed 
galaxies, but for Ab < 0, the average bias is closer to unity, 
and the prior is better matched. As a consequence, it is not 
surprising that the modelling leads to quite large biases in 
the parameter estimates for Ab > 0. 

In the analysis so far, we have assumed that we know 
the slope of the number counts exactly. In Fig. [6] we show 
how sensitive the errors, both statistical and systematic, are 
to errors in the number count slope. We assume there is 
a gaussian random scatter with width cr^ about the true 
value of a in each redshift bin. Note that these panels have 



one logarithmic axis. We see that if we null and assume that 
there is no residual k6 term (the 'naive nulling' curves), then 
scatter around the fiducial number count slope is very seri- 
ous for parameter errors, if the scatter exceeds about 0.02. 
For nulling where we also model the contaminating k6 term, 
the situation is much better, with parameter estimation after 
marginalisation over the galaxy bias being rather insensitive 
to variations of the slope. In the worst cases, the parameter 
bias is of order the systematic error, and if the slope can be 
constrained to an error of less than 0.02, the systematic er- 
ror is subdominant for all parameters. Note that the rather 
strange behaviour of the bias in the spectral index Us for the 
nulled case is due to random fluctuations from a relatively 
small number of trials, but it is in any case subdominant for 
all values of the scatter probed. 
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Figure 5. Statistical la errors (black lines) and biases (grey /red lines) for different cosmological parameters as a function of the offset 
Ab of the galaxy bias model, keeping all other parameters at their fiducial values. Solid curves correspond to the results for standard 
analysis including marginalisation over the galaxy bias, dotted curves to the combined nulling and galaxy bias marginalisation approach. 



5 DISCUSSION 

Observed galaxy clustering for deep surveys is determined by 
a combination of intrinsic clustering and changes in number 
density due to the magnification and amplification effects 
of galaxy lensing. If the lensing effect could be isolated, it 
could be a useful cosmological probe, as it is independent 
of galaxy bias, and depends rather cleanly on cosmological 
parameters. The difficulty is that the signal is small in com- 
parison with the intrinsic clustering term. In this paper, we 
have presented a method to isolate as far as possible the part 
of the clustering signal which is due to cosmic magnification. 
This involves exploiting the known redshift dependence of 
the magnification signal to remove the cross-terms which 
link the magnification with foreground intrinsic clustering. 
By isolating the magnification term, we aim to remove the 
dependence of clustering on the galaxy bias, which is subject 
to uncertainties connected with galaxy formation efficiency. 



leaving a clean test of cosmology relying only on the be- 
haviour of the matter power spectrum and distance-redshift 
relations, both of which are much more easily predictable 
theoretically than galaxy bias. This is a very challenging 
task, because the contamination of the clean signal by intrin- 
sic clustering is very large, and we have been partially suc- 
cessful. The intrinsic clustering signal can be reduced very 
substantially by the 'nulling' technique we have employed 
(by a factor of about 100), but not to a negligible level, so 
the residual contamination still has to be modelled. How- 
ever, since the amplitude of the residual contamination is 
much smaller than the signal before nulling, we do not need 
to model bias as accurately. By investigating an example 
cosmological setup where galaxy bias has a mild evolution 
with redshift, and modelling it incorrectly as a constant bias, 
we find that full modelling of the clustering signal (includ- 
ing the intrinsic clustering term) has small statistical errors. 



10 Alan F. Heavens and Benjamin Joachimi 



10 



0.1 



0.01 



Q 






a; naive nulling 




b; naive nulling 




a; nulled 




b; nulled 





0.02 0.04 



0.06 0.08 



10 



0.1 



0.01 



a; naive nulling 
b; naive nulling 
a; nulled 
b; nulled 



0.1 0.12 0.14 



0.02 0.04 0.06 0.08 0.1 0.12 0.14 



10 



0.1 



0.01 



a; naive nulling 
b; naive nulling 
a; nulled 
b; nulled 



0.02 0.04 



0.06 0.08 



10 



0.1 



0.01 



■ ^0 










a; naive nulling 




b; naive nulling 




a; nulled 




b; nulled 







0.1 0.12 0.14 



0.02 0.04 



0.06 0.08 



0.1 0.12 0.14 



Figure 6. Statistical la errors (black lines) and biases (red/grey lines) for different cosmological parameters as a function of the la 
scatter aa of the values of a in the different redshift bins, keeping all other parameters at their fiducial values. Chain-dotted curves 
correspond to the results for naive nulling (without marginalisation over the galaxy bias), dotted curves to the combined nulling and 
galaxy bias marginalisation approach. The values of a are drawn from a Gaussian distribution centred on the fiducial values of a and 
with a scatter a^- The results plotted are averaged over 10 realisations. To avoid noise, the same seed is used for each value of a^- 



principally because the kS term is carrying quite a lot of cos- 
mological information, but the incorrect parametrisation of 
the galaxy bias leads to systematic errors in the cosmologi- 
cal parameters which are typically significantly larger than 
the statistical errors. Nulling, on the other hand, reduces the 
sensitivity to the galaxy bias modelling, as most of the effect 
of galaxy bias is removed. The result is that the systematic 
errors in the cosmological parameter estimates are reduced 
to a low level, at the expense of increases in statistical er- 
rors, which in some cases are quite large. The reason for the 
increase in error bars is that the term is relatively large 
(in contrast to the GI term in cosmic shear), and it contains 
cosmological information, so removing it does increase the 
errors. Thus nulling gives a robust alternative analysis to 
modelling the full clustering signal, as it is less sensitive to 
modelling of the galaxy bias, but it is a conservative anal- 



ysis. In short, full modelling is subject to a bias of a size 
which we may not know, as it depends on how wrong our 
modelling of galaxy bias is. Nulling protects us against this. 
However, the realistic position is not so pessimistic, as we do 
in principle have empirical measurements of bias available, 
from the auto-correlation of number counts in tomographic 
bins, or from the galaxy bispectrum, which may prove help- 
ful. 

One interesting technical point is how we are able to im- 
prove our parameter estimation by throwing away data. We 
see that the weighted cross-powers are simply linear combi- 
nations of individual bin-bin cross-power, but we have sim- 
ply removed one linear combination of the bin number den- 
sities. The reason is as follows. The removed combination 
has a very large contribution from intrinsic clustering, and 
is thus sensitive to how the galaxy bias is modelled. If this 
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is done imperfectly, for example by assuming some redshift 
dependence which does not hold in nature, then the pa- 
rameter estimation process may try to fit this combination, 
at the expense of a poorer fit to the nulled combinations. 
As a result the parameter estimation may be systematically 
wrong. By nulling the intrinsic signal, the cross-power spec- 
tra have a much greater signal coming from magnification, 
independent of bias, so the fitting is more likely to find the 
correct parameters. Note that if we know the bias evolution 
precisely, then it will always be better to model all the cross- 
power, as this does not lose information. In practical cases, 
it is an open question which method is preferred, and this 
will be the subject of a further paper. 

The key issue for the use of magnification for cosmo- 
logical parameter estimation will be the level of systematic 
errors, since surveys which are beginning now and which are 
planned for the future are large enough that it is likely that 
the statistical errors will become irrelevant. We have not ad- 
dressed this issue in the paper, except to an extent discussing 
the systematic errors arising from imperfect bias modelling. 
Evidently, zero-point magnitude errors will feed into the er- 
ror budget, although recent improvements in photometry 
(e.g. Ilbert, O., et al.|[2006| are encouraging. Note that the 



principal effect here will be angular on the sky, so the red- 
shift dependence of the lensing signal may help to alleviate 
this issue. Perhaps of most concern is the photometric red- 
shift estimation (e.g. Hildebrandt et al._2010 and references 
therein), where catastrophic outliers will correlate intrin- 
sic clustering over large redshift separations, and mimic the 
magnification signal to some degree, and the nulling process 



itself is also sensitive to outliers (see Joachimi Schnei- 
[der||2QQ 9 for details). Both of these require further study to 
see what requirements will be needed to be met for cosmic 
magnification to be a viable technique. 
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